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V ' Abstract 
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. The data recorded at a centre-of-mass energy of 189 GeV by the OPAL detector at LEP are used to 
^ [ search for trihnear couplings of the neutral gauge bosons in the process e'''e~— > Z7. The cross-sections 
for multihadronic events with an energetic isolated photon, and for events with a high energy photon 
accompanied by missing energy are measured. These cross-sections and the photon energy, polar 
angle and isolation angle distributions are compared to the Standard Model predictions and to the 
theoretical expectations of a model allowing for Z7Z and Z77 vertices. Since no significant deviations 
with respect to the Standard Model expectations are found, constraints are derived on the strength 
of neutral trilinear gauge couplings. 
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1 Introduction 



The self-interactions of the gauge bosons are consequences of the non-Abehan structure of the elec- 
troweak sector of the Standard Model; therefore, the strength of trilinear and quartic gauge couplings 
is predicted as a result of the gauge symmetry of the theory. The study of trilinear gauge boson cou- 
plings in two-boson production processes is within the reach of existing accelerators and measurements 
of WWZ and WW7 couplings are being performed with increasing precision in e'''e~ j^] and pp ||^] 
collisions. While non-zero values are predicted for the couplings in the triple and quartic vertices 
involving charged gauge bosons, the tree level vertices Z7Z, Z77 and ZZZ are not generated by the 
Standard Model Lagrangian; higher order corrections through virtual loops contribute at the level of 
10~^ 1^, well below the current experimental sensitivity. Nevertheless, new phenomena with a charac- 
teristic mass scale above the present experimental threshold might lead to tree-level neutral trilinear 
gauge couplings (NTGC) in the effective Lagrangian parametrising the residual low energy effects 
from new physics. For example, as suggested in |^], virtual effects from new heavy fermions having 
non-standard couplings to the gauge bosons might generate sizeable anomalous contributions. 

The most general Z7V vertex (where V is the intermediate virtual boson, either photon or a Z) 
compatible with Lorentz invariance and electro-magnetic gauge invariance involves four independent 
operators, corresponding to the allowed helicity states for the Z7 pair Therefore, in a model 

independent description, there exist eight couplings: four of them {hf, i = 1, •..,4) corresponding to 
V=Z and four (hj ) corresponding to V=7. The vertex function was first given in [Q. In this analysis, 
we adopt the most recent convention established in The lowest dimensional operators associated 
to the h^''^ {i = 1,3) couplings are of dimension six, while dimension eight operators are associated to 
h^'^ {i = 2,4). As discussed in Q, no additional symmetry constraints, such as the SU(2)xU(l) gauge 
invariance usually assumed in the case of WWV anomalous couplings j^, can help in reducing the 
number of free parameters in the neutral gauge vertex, since operators of even higher dimensionality 
would be required. Therefore a model independent approach is used which retains all eight couplings. 

In this paper, the process e+e~— > Z7 at ^/s = 189 GeV is investigated through the final states 
qq7 and i^i^'j, the dominant decay modes of the Z boson, with the aim of searching for Z7Z and 
Z77 couplings. Experimental constraints on these couplings have been produced in the past from the 



analysis of LEP data at lower energies |1C] and of the Tevatron data [11|. The analysis presented 
here has higher sensitivity due to the increased centre-of-mass energy and due to the large data 
sample collected during the 1998 operation of LEP. It also benefits from the recent clarification |^] 
of the theoretical framework in which neutral gauge boson self-interactions can be described. A 



recent analysis based on data collected at ^/s = 189 GeV by the L3 collaboration [12| adopts the same 
convention as in this paper. On the other hand, due to this different convention, the comparison of 
the results presented here and in [12| with previous published results is not straightforward. 



In e^e^ collisions, the production of Z7 final states via anomalous neutral gauge couplings has 
a large irreducible Standard Model background from production with hard initial state radiation 
(ISR). Small contributions to qq7 production also arise from e''~e~— > 77* — > 7qq and from final state 
radiation in qq production. In the I'V^ channel, final states arising from the exchange of a W boson in 
the t-channel also contribute to the cross-section, although their contribution is small within the Z7 
signal acceptance used in this analysis. As a general property, the Z and 7 produced at the anomalous 
vertices are more isotropically distributed than the dominant Standard Model background, which is 
characterized by the strongly forward peaked angular distribution of initial state radiation. Therefore, 
deviations from the Standard Model predictions due to Z7V couplings would be more pronounced 
for large angles between the beam direction and the photon. In the qq7 final state, the angular 
distribution of the jets can also be exploited in order to gain sensitivity to anomalous couplings, due 
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to the resulting enhancement of the longitudinal polarisation of the Z boson affecting the fermion 
decay angle. On the other hand, the photon energy spectrum has a marginal sensitivity, due to the 
kinematic constraints from the fixed centre-of-mass energy and the narrow Z resonance. Since all the 
terms in the Z7V vertex are proportional to the momenta of the gauge bosons involved, this results 
in an enhancement of the sensitivity to NTGC as the centre-of-mass energy increases. Finally, the 
experimental signature of the anomalous h^'^ couplings depends on the CP parity of the associated 
operators. The vertex terms proportional to h^'^ and h^'^ violate CP and, hence, do not interfere 
with the CP conserving Standard Model amplitudes. As a result the total and differential cross- 
sections receive only additive contributions from the anomalous processes. The remaining couplings, 
associated to CP even terms, lead to amplitudes which interfere with the Standard Model; therefore 
the differential and total cross-sections are enhanced or suppressed depending on the sign and the size 
of the h^'^ {i = 3, 4) couplings. 



2 Detector and Monte Carlo Simulation 



The OPAL detector, described in detail in |13], consists of a central tracking system inside a solenoid 
providing a magnetic field of 0.435 T, and of an electromagnetic calorimeter, complemented by a 
presampling system and an array of scintillation counters for time-of- flight measurements; hadron 
calorimetry is obtained by instrumenting the magnet return yoke which is surrounded by muon cham- 
bers. A system of forward calorimeters extends the angular coverage of the detector down to a polar 
angl^ of 24 mrad. However, due to the installation in 1996 of a thick tungsten shield designed to 
protect the tracking chambers from synchrotron radiation background, the effective limit of electro- 
magnetic hermeticity is around 33 mrad. The integrated luminosity of the data samples is determined 
from the rate of small angle Bhabha scattering events observed in the silicon-tungsten calorimeter [0] 
with a precision of 0.22%. 



Track reconstruction is performed by combining the information from a silicon microvertex detec- 
tor, a vertex drift chamber, a large volume jet drift chamber and an outer layer of drift chambers 
for the measurement of the z coordinate. The most relevant subdetector for the event topologies 
used in the analysis presented here is the electromagnetic calorimeter. It consists of an array of 9440 
lead-glass blocks in the barrel (| cos 9\ < 0.82) arranged in an almost-pointing geometry and two dome- 
shaped end caps, each of 1132 longitudinally aligned lead-glass blocks, covering the polar angle range 
0.81 < I cos^l < 0.984. Trigger signals |l^, based on energy deposits in the lead-glass blocks and also 
on a coincidence of energy in the barrel electromagnetic calorimeter and a hit in the time-of-flight 
system, guarantee full trigger efficiency for both the qq7 and the uV'j events falling within the signal 
definition criteria used in this analysis. 



The Standard Model processes leading to the qq7 and uV^y final states have been simulated using 
the Monte Carlo generators KK2f [|l6| for the hadronic channel and KORALZ |l^ and Nunugpv98 @ 
for the missing energy channel. For the qq7 channel, the fragmentation, which includes photon radia- 
tion from the quarks, and the hadronization are simulated with the Jetset package [19| tuned on the 
basis of extensive studies of hadronic events at the Z resonance as described in |2^. The grc4f [^ ] 
generator has been used to estimate the background to the qq7 channel from four-fermion produc- 
tion. The contribution to the background from two-photon interactions has been studied with the 



Monte Carlo generators Phojet [22|, for the untagged and double-tagged events, and Herwig 
for the single-tagged events and charged current deep inelastic scattering events. In the i^Vj channel, 

^In the OPAL coordinate system, 9 is the polar angle defined with respect to the electron beam direction and is 
the azimuthal angle. 
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the contamination from Bhabha events has been estimated using the Bhwide |24| and Teegg 



generators. The contamination from four-fermion production has been studied using grc4f and 



KORALW 1 26], while to determine the background from two-photon interactions the Vermaseren [^] 
Monte Carlo generator has been used. The Radcor Monte Carlo has been used to study the 
e+e~ 77 background and the energy response of the calorimeter to photons. All the Monte Carlo 
samples described above were processed through the OPAL detector simulation |29]. For the inter- 
pretation of the data, as will be described in more detail in section ^, a Monte Carlo generator 
based on the matrix element for ff7 production in e'^e" collisions and including the contributions 
from NTGC, has been used. 



3 Event Selection and Cross-Section Measurements 
3.1 The selection of qq7 events 

The selection of hadronic events with isolated high energy photons is performed on events preselected 
as high multiplicity hadronic events in a data sample corresponding to an integrated luminosity of 176.2 
pb~^ with an average centre-of-mass energy of 188.6 GeV. The preselection criteria, described in 
are based on the track and cluster multiplicity, on the visible energy and on the longitudinal imbalance 
of the energy measured in the electromagnetic calorimeter. The events satisfying the preselection 
requirements are processed by a photon search algorithm. 

The photon identification is based on an algorithm optimised for photon search in hadronic events 



described in [32|. Electromagnetic clusters without associated tracks in the central detector are ac- 
cepted as photon candidates if their energy is higher than 5% of the beam energy and their polar angle 
lies in the acceptance region of the lead-glass calorimeter. The number of lead-glass blocks involved 
and the energy sharing among them are required to correspond to typical patterns defined for photon 
identification in the OPAL calorimeter. An isolation criterion is then applied in order to reject electro- 
magnetic clusters associated with jets. The total energy deposition in the electromagnetic calorimeter 
(not associated to the photon candidate) within a cone of 15° around the photon flight direction is 
required to be less than 2 GeV. In addition, the sum of the momenta of tracks which, extrapolated 
to the calorimeter surface, fall inside a 15° cone around the photon impact point is also required to 
be lower than 2 GeV. Finally, systems of one or two well reconstructed tracks associated with elec- 
tromagnetic clusters, consistent with a photon conversion according to the criteria described in 
are included if the reconstructed photon satisfies the criteria listed above. The photon identification 
algorithm has been extensively studied in order to assess the level of accuracy of the estimate of the 
efficiency obtained from Monte Carlo qq7 events. In particular, the selection criteria of the algorithm 
have been adjusted in order to minimise their sensitivity to unsatisfactory modelling in the simulation. 
A residual discrepancy in the identification efficiency for converted photons has been observed and 
taken into account as a correction of 1.25% to the overall selection efficiency. As an example of the 
quality of the modelling of the photon identification algorithm, figure |^ shows the angle between the 
photon candidate and the closest track in the event and the total charged energy in the isolation cone. 

After the photon search, all the clusters and tracks in the event which are not associated to the 



most energetic photon candidate are grouped into jets according to the Durham fcy |34] scheme with 
resolution parameter y = 0.02. If more than four jets are reconstructed, the event is forced to have 
four jets in addition to the isolated high energy photon. 

In e"'"e~— > Z7 at ^/s = 189 GeV the photon energy spectrum is peaked at approximately 72 GeV, 
reflecting the sharp Z resonance. In order to select the topology corresponding to a high energy photon 
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recoiling against a hadronic system of invariant mass equal to the Z boson mass, the signal definition 
is based on kinematic cuts applied to the most energetic photon in the event: 

50 GeV < < 90 GeV ; 
15° <9^< 165°, 

where and 0^ are the photon energy and polar angle, respectively. 

To improve the photon energy resolution and suppress further the background from non-qq7 events 
and the feedthrough from events outside the signal definition, a kinematic fit is applied to all the events 
with at least one photon of energy larger than 30 GeV. The fit imposes energy and momentum 
conservation using as input the photon and the jet momenta. Undetected ISR is allowed to compensate 
for missing longitudinal momentum in the beam pipe region if the probability of the fit is smaller 
than 1%. 

From a study of Monte Carlo qq7 events, the kinematic fit improves the photon energy resolution by 
a factor two. The events for which the fit converges (99.5%) are finally selected if the fitted values of 
and 6^ satisfy the signal definition cuts. In order to suppress contamination from photons originating 
from the jets, only the sub-sample in which a-y-jct > 30°, where a^_jet is the angle between the photon 
and the closest jet, is retained for the analysis. 

The selection efficiency and the feedthrough in the kinematic acceptance are listed in table |. They 
have been estimated from fully simulated e"'"e~— >qq Monte Carlo events, where the signal-like topology 
arises from radiative return to the Z resonance. The efficiency has been calculated with respect to the 
kinematic signal acceptance defined above. The feedthrough represents the fraction of selected events 
which do not belong to the kinematic acceptance; it is mainly due to resolution effects, but it includes 
also a small contamination (0.35%) due to non-ISR photon candidates in e^e~— >qq events falling in 
the signal acceptance after reconstruction. Both the efficiency and the feedthrough are corrected for 
the aforementioned residual disagreement between the performance of the photon search algorithm in 
the data and in the Monte Carlo. 

The numbers of events selected in the data and in the Monte Carlo samples are listed in table H. 
Figure ^ shows the distributions of the photon energy and polar angle and angular separation with 
respect to the closest jet for the events selected in the data, compared with the Standard Model 
expectation from Monte Carlo. The total background (2.5% of the selected events) comes from four- 
fermion production (1.59%) two-photon interactions (0.59%) and tau-pair production (0.34%). The 
agreement between the data and the Monte Carlo is in general satisfactory, except perhaps in the 
photon energy distribution where a slight deficit of events is observed in the radiative return peak. 

3.2 The selection of wV^ events 

The selection of events with an isolated high energy photon accompanied by missing energy and low 
activity in the detector follows the single-photon analysis described in ||3^. The data sample used in 
the analysis corresponds to an integrated luminosity of 177.3 pb~^, with an average centre-of-mass 
energy of 188.6 GeV. After the single photon selection, 643 events are retained in the data. The 
same additional conditions which define the kinematic acceptance in the qq7 channel, 50 GeV < E^ < 
90 GeV and 15° < 6^ < 165°, are then applied to the most energetic photon in the electromagnetic 
calorimeter. 
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The predicted efficiency and the feedthrough in the kinematic acceptance are hsted in table ||. 
They have been estimated averaging the predictions of the KORALZ and the Nunugpv98 Monte 
Carlo, which in the signal acceptance agree within (0.2ib0.5)% in the efficiency and (lib9)% in the 
feedthrough. 

The number of events selected in the data are listed in table |^ together with those expected from 
all the relevant physics processes and from instrumental backgrounds. As potential sources of physics 
background, four-fermion processes, radiative Bhabha and e~''e~ — > 77 events have been considered, 
while the residual cosmic ray and beam-related contamination have been estimated using control 
samples enriched in these backgrounds. Both these sources result in an overall negligible (0.24%) 
contribution to the selected events. Figure |^ shows the distributions of the energy and of the polar 
angle of the photon for the selected events, compared with the Standard Model expectation. 

3.3 Cross-section Measurements 

Prom the number of observed events in each channel and from the predicted efficiency, feedthrough 
and backgrounds, as presented in sections ^1 and |3[2, the cross-sections within the kinematic signal 
acceptance are measured to be: 

o-qq^ = 9.42 ± 0.25 (stat.) ±0.15 (syst.) pb 
a^u^ = 2.52 ±0.13 (stat.) ±0.05 (syst.) pb. 

These measurements are in reasonable agreement, within the errors which are dominated by the 
statistical uncertainty, with the predictions from the Standard Model, which are respectively Cq^ = 
9.75 ± 0.03 pb and a^f^ = 2.81 ± 0.02 pb. The predicted qq7 cross-section is obtained from KK2f, 
while the uv^ cross-section is the average of the predictions of KORALZ and NuNUGPV98 which are 
consistent within (1.6 ± 1.1)%. The errors associated with the Standard Model predictions come from 
the Monte Carlo statistics. Taking into account the experimental systematic uncertainties discussed 
in the following section, the overall discrepancy corresponds to —1.1 standard deviations for the 
qq7 channel and —2.1 for the vV"f channel. 

3.4 Systematic Errors 

The different sources of systematic uncertainties affecting the cross-section measurements are sum- 
marised in table |^ and are discussed in the following: 

Systematics specific to the qq7 channel: 

- The uncertainty on the selection efficiency. Contributions to this uncertainty come from imper- 
fect modelling of the material in the detector, affecting the estimate of the photon conversion rate 
(0.77%), from inaccuracies in the modelling of the identification efficiency for converted photons 
(0.62%), from uncertainty in the simulation of the track-cluster association at forward angles 
(0.64%) and from the photon isolation criterion (0.26%). The total uncertainty coming from the 
modelling of the selection efficiency translates into a relative error on the cross-section of 1.2%. 
The uncertainties arising from limited Monte Carlo statistics used to evaluate the efficiency and 
the feedthrough are respectively 0.14% and 0.07%. 
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- The modelling of the jet reconstruction. An additional smearing of the jet energies and directions 
and a shift of the jet energy scale has been applied in the Monte Carlo, on the basis of an extensive 
comparison of two-jets events in calibration data collected at the peak and in the simulation. 
The resulting variation (0.36%) in the cross-section has been assigned as a systematic error. 

- The sensitivity of the analysis to the jet multiplicity in the event. The results of a modified 
analysis, where each event has been forced to contain exactly two jets in addition to the isolated 
photon, have been compared with those from the standard analysis. The difference has been 
assigned as a systematic error (0.69%). 

- The uncertainty related to the models used to simulate the hadronisation process. This uncer- 
tainty (0.08%) has been evaluated comparing the results when either the Jetset or the Herwig 
hadronisation schemes have been used in the Monte Carlo. 

- The uncertainty (0.61%) due to the background subtraction; this is dominated by the 100% 
uncertainty on the normalisation of the background from two-photon interactions as predicted 
by Herwig, which is expected to give the best description of two-photon interactions in the 
data, and by Pythia. This large uncertainty is assigned to cover possible mismodelling of the 
hard-fragmentation processes in the very small fraction of the two-photon cross-section retained 
in the selection, as suggested by a comparison with the F2GEN |^6| generator. Finally, the 
uncertainty related to the contamination from fake ISR photons in non-radiative qq events is 
estimated to be 0.35%. 



Systematics specific to the vv^ channel: 



- The uncertainty on the selection efficiency. Systematic errors come from the estimate of the 
detector occupancy (1%) and from the imperfect description of the conversion probability and 
of the conversion tagging efficiency in the Monte Carlo (0.8%). A 0.9% uncertainty is assigned 
to account for systematic effects in the simulation of the cuts used to reject radiative Bhabha 
and e''"e~ — > 77 events. The uncertainty on the efficiency of the timing cuts used to reject 
background from cosmic rays translates into a systematic error of 0.5%, while an uncertainty 
of 0.6% is assigned to cover possible mismodelling of the other cuts complementing the timing 
requirements in the cosmic and instrumental background rejection. The overall uncertainty 
on the modelling of the selection efficiency is estimated to be 1.7%. Finally, the statistical 
uncertainty on the efficiency and feedthrough as estimated from Monte Carlo contributes with 
a 0.3% uncertainty on the cross-section. 

- The uncertainty on the background from processes other than vV^ , amounting to 0.25%. This 
includes a 50% systematic uncertainty on the expected background, which was assigned to cover 
possible mismodelling of the vetoes used in the rejection. 



Systematic uncertainties common to both analyses: 



The uncertainty on the angular acceptance arising from residual biasses in the coordinate recon- 
struction and from the absolute knowledge of the detector geometry. The overall uncertainty 
on the position of the electromagnetic showers, estimated in |37| to be 0.001 rad, results in a 
relative uncertainty of 0.19%. 

The effects of mismodelling of the energy scale and the resolution of the electromagnetic calorime- 
ter. They have been investigated using a control sample of about one thousand strongly collinear 
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e"'"e~ — > 77 events, whose energies are expected to be very close to the beam energy. The Monte 
Carlo has been found to reproduce the electromagnetic calorimeter energy scale and resolution 
in the data within respectively 0.3% and 10%. These differences are in very good agreement 
with those observed using as a reference process a sample of wide angle Bhabha events |37|. The 



systematic error has then been estimated by modifying the absolute energy scale and resolution 
in the Monte Carlo within these limits and the corresponding variation of the cross-section has 
been taken as systematic error. 

The uncertainty on the beam energy (±20 MeV). The effect on the cross-section was evaluated 
by appropriately scaling the energy of the most energetic photon in the Monte Carlo, while 
leaving unchanged the invariant mass of the system recoiling against the photon. 

The uncertainty on the measurement of the integrated luminosity, 0.22%. 



4 Data Interpretation 

Since the data presented in the previous section show no evidence of deviations with respect to the 
Standard Model expectation, they are used to derive bounds on the strength of anomalous Z7Z and 
Z77 couplings. 



4.1 Analysis Procedure 



The analysis is based on a comparison of the measured event rate and of the differential distributions 
with the theoretical predictions for the Standard Model processes and possible contributions from 
NTGC. The energy spectrum and the cosine of the polar angle of the photon are used in both channels, 
while in the qq7 channel the cosine of the angle between the photon and the closest jet is also used. 
Due to the almost monochromatic photon energy spectrum, this angle is strongly correlated to the 
quark emission angle in the Z decay rest frame, which is sensitive to NTGC. 

The theoretical predictions as a function of different values of the anomalous couplings h^''*' are 
obtained from a Monte Carlo generator |3C] for ff7 production in e+e~ collisions. This generator is 
based on the full matrix element, in the lowest order approximation, for all the relevant Standard 
Model processes and for processes generated by anomalous trilinear neutral gauge couplings. The 
only Standard Model contribution missing in the calculation is the t-channel W boson exchange in the 
vV^ channel, which has been estimated from NuNUGPv98 and included in the reweighting procedure 
described in the following. The effect of higher order QED corrections from initial state radiation, 
which have been found to reduce the contribution from anomalous couplings by typically 15%, has 
been incorporated into the calculation using a collinear radiator function from the Excalibur |^] 
Monte Carlo. 



The events selected in the data have been classified in 5 x 4 x 4 unequal bins of the three-dimensional 
(E^, cos 6^, cos a^_jet) space. In the case of the vV"f channel, the angle between the photon and the Z 
decay products is not experimentally accessible and therefore integrated out. The NTGC-dependent 
theoretical prediction for the population in each cell, provided by the Monte Carlo calculation, has 
been modified to allow for the reconstruction efficiency and resolution effects, as determined from 
a large sample of fully simulated Z''/7 qq and 71^ — > vv Standard Model Monte Carlo events. 
The total number of expected events and the population of each cell as a function of the anomalous 
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couplings, N(hj), is determined by reweighting the number Nsm of accepted events predicted by the 
fully simulated Standard Model Monte Carlo according to: 

N(h,)=NsM(l + 5(h.)), ^(h,) = N^^c(h^) - N.ec(h^^ ^ 

where Nrec is the number of reconstructed events from the NTGC-dependent theoretical prediction, 
modified for efficiency and resolution as explained above. Figures ^ and ^ show how the distributions of 
the kinematic variables, folded with detector effects, are modified by a particular choice of anomalous 
couplings in the qq7 and vV^ channels. In order to disentangle the effect of NTGC on the event rate 
and on the differential distributions the latter have been normalised to the number of events selected 
in the data. 

The theoretical expectations for the event rate and the differential distributions have been fitted 
to the data independently for the two channels, under the hypothesis that only one coupling at a time 
is non-zero. The most probable values of the anomalous couplings are determined by minimising the 
negative log-likelihood defined as follows: 

-LogL = -LogP(N°b^N(h)) - ^LogP(Nf ^N,(h)) 

i 

where P(N°'^**, N(h)) is the Poisson probability of observing the number of events N°''** if the expectation 
is N(h). The index j runs over the number of cells defining the multidimensional distributions and the 
condition ^jNj(h) = N"*^*^ is imposed to disentangle the contributions to the likelihood from event 
rate and distributions. 

The fit procedure has been tested on Monte Carlo samples of Standard Model qq7 and uv^ 
events reweighted for NTGC effects. The central values of the fit results have been found to correctly 
reproduce the values of the input NTGC parameters. In order to check the reliability of the errors 
on the fit results, tests have been performed on several Standard Model Monte Carlo samples of 
size corresponding to the data luminosity. The distributions of the central values of the couplings 
determined by the fit are found to be consistent with those expected from the statistical sensitivity. 



4.2 Results on Trilinear Neutral Gauge Couplings 

The values of the anomalous couplings and the statistical errors obtained from the likelihood fit of the 
qq7 and vu^ event rate and distributions are listed in table |5[ together with the expected statistical 
sensitivity of the analysis which would be achieved in the case of perfect agreement between data 
and the Standard Model predictions. In general, the total event rate and the differential distributions 
have similar sensitivities to NTGC. However, due to the quadratic dependence of the cross-section 
on the couplings, the fit of the event rate can only determine the value of the couplings with a two- 
fold ambiguity. In the case of the CP violating couplings, h^'"' {i = 1,2), which lead to amplitudes 
that do not interfere with the Standard Model amplitudes, the sign of the couplings is completely 
undefined, but both the cross-section and the distributions provide a determination of the absolute 
value of the couplings. In the case of the CP conserving couplings, h^'^ (i = 3,4), the interference 
with the Standard Model amplitudes results in distributions of the kinematic variables which produce 
a unique minimum in the associated — log L function, thus removing the two-fold ambiguity arising 
from the event rate information. 

The 95% Confidence Level (C.L.) bounds on the eight anomalous couplings have been obtained 
by convolving the likelihood function with a Gaussian whose width a corresponds to the systematic 
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uncertainty on the individual parameters, as estimated in section The central values and the 
95% C.L. intervals resulting from the combination of the two channels are given in table ^. Figures 
^ and |7| show the corresponding negative log-likelihood curves for the individual channels and their 
combination. In combining the results, the correlations between the systematic uncertainties have 
been taken into account. 



4.3 Systematic Errors 



The impact of several sources of systematic uncertainty on the NTGC determination has been assessed. 
Most of the effects considered have already been discussed in the context of the systematic error on 
the cross-section measurements. A few more are related to the modelling of NTGC effects on the 
observables used in the likelihood fit and to the reference Standard Model predictions. For all the 
sources of systematic uncertainty, a symmetric error is assigned based on the maximum absolute shift 
between the central value, the lower and the upper edges of the 68% C.L. interval as obtained in the 
standard analysis (table ^) and in a specific fit to the data performed to simulate the systematic effect. 
The dominant sources of systematic uncertainty come from: 



- The effects of the modelling of the selection efficiency. These have been evaluated using the same 



methods as discussed in section p.4| , and have been assigned as scale uncertainties on the expected 
number of events in each channel. Where relevant, the angular dependence of the uncertainties 
has been taken into account. For the qq7 channel, this includes also the uncertainties coming 
from the jet multiplicity, the hadronisation modelling and the jet parameter smearing. 

The theoretical uncertainty on the Standard Model prediction. Considerations about missing 
higher order corrections in the KK2f Monte Carlo lead to an estimate of the theoretical uncer- 
tainty of the order of 1% in the qq7 channel | If: ] . In the vV^ channel, a 2% theoretical uncertainty 



has been assigned; this uncertainty covers the observed differences between the kinematic cross- 
sections and the selection efficiencies estimated by KORALZ and NuNUGPv98, and is consistent 
with a recent comparison between different theoretical calculations presented in [^. These 
theoretical uncertainties have been assigned as an overall normalisation error on the Standard 
Model prediction in each channel. 



Other minor contributions to the systematic error have been considered: 



- The contribution of the background to the fit result. This has been separated into a component 
due to a scale factor (corresponding to the Monte Carlo statistical uncertainty on the background 
absolute rate) and a component due to the modelling of the background shape. The error due 
to the background normalisation is evaluated as the maximum effect observed when increasing 
or decreasing the total background rate by one standard deviation. The error due to the shape 
modelling has been conservatively assessed by assuming a fiat distribution of the background 
events in the signal acceptance region. In the vV^y channel, due to the very small contamination, 
only the uncertainty on the background rate has been considered. 

- The uncertainty related to the Monte Carlo statistics. This effect has been estimated by applying 
the likelihood fit to the data using, as Standard Model prediction in each bin, a number of events 
generated according to a Poisson distribution with average equal to the prediction of the reference 
Monte Carlo sample and rescaled to the integrated luminosity of the data. The maximum among 
the r.m.s of the central values and of the lower and upper 68% C.L. limits on each coupling, over 
a thousand fits, has been assigned as systematic error. 
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The uncertainty on the reweighting procedure coming from the hmited Monte Carlo statistics 
used to evaluate the corrections for resolution and selection efficiency, and from the limited 
Monte Carlo statistics used in the calculation of the generator- level weights. 

The uncertainty on the missing t-channel W exchange contribution, which applies only to the 
vV^ channel. The procedure used to correct for the missing W t-channel in the reweighting 
procedure accounts for this process only in the pure Standard Model contribution, but not in 
the interference between the Standard Model process and the process leading to anomalous 
coupling. Since the W t-channel is calculated to contribute less than 4% to the Standard Model 
cross-section in the signal acceptance, the effect of it being neglected in the interference is 
expected to be of order 2%. The interference term has been conservatively varied by ±4% and 
the differences in the fit results have been assigned as a systematic error. 

The systematic errors arising from the calorimeter energy scale and resolution, the modelling of 
the angular cut, the uncertainty on the beam energy and on the integrated luminosity. They 



have been evaluated with the same method as discussed in section 3.4 and have been treated as 
fully correlated between the qq7 channel and the vV^ channel. 



The systematic errors on the couplings are summarised in table |^ for both the qq7 channel and 
the vv^ channel. For all the couplings the total systematic error is small compared to the statistical 
uncertainty, except for hg^, where the size of the systematic error reaches 60% of the statistical error. 



5 Conclusions 



Using the data collected at \fs = 189 GeV, the cross-sections and the differential distributions for 
hadronic events with a high energy isolated photon and for events with an energetic photon and 
missing energy have been measured to search for possible contributions from anomalous Z7Z and 
Z77 couplings. Since no significant evidence of deviations with respect to the Standard Model is 
observed, 95% C.L. limits on the eight trilinear neutral gauge couplings h^'^ have been derived. These 
limits do not yet allow to place significant constraints on specific models of new physics leading to 
effective anomalous couplings in the neutral sector. Nevertheless, these results and those presented 
in , which are compatible and of equivalent sensitivity, represent the best available investigations in 
e''~e~ — > Z7 of the neutral gauge boson self- interactions, in the recently revised theoretical framework 
describing the general Z7V (V = Z,7) vertex. 
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Efficiencies and Feedtiirougli 


Cfiannel 


EfRciency (%) 


Feedtfirougfi (%) 


qq7 


87.58 ±0.12 
81.04 ±0.22 


2.21 ±0.07 
1.80 ±0.08 



Table 1: Tfie efficiency and tfie feedtfirougfi for tfie qq7 and vV^ selections. In the uV^ channel, they 
are the averages of the KORALZ and the NUNUGPV Monte Carlo. The errors are from Monte Carlo 
statistics. 



Selected Events 


Observed 


qq7 SM 


Four-Fermion 


Two-Photon 




Total expected 


1525 


1538.7 ±5.5 


24.3 ± 0.9 


9.0 ± 1.0 


5.2 ±0.3 


1577.2 ± 5.7 



Table 2: The number of selected qq7 events in the data, compared with the Monte Carlo expectations 
for the signal and for the background, normalised to the data integrated luminosity. The errors are 
from Monte Carlo statistics. 



Selected Events 


Observed 


uV^ SM 


Physics Background 


Instrumental Background 


Total expected 


370 


411.6 ±2.5 


0.72 ± 0.83 


0.18 ±0.18 


412.5 ±2.6 



Table 3: The number of vV^ events selected in data, the expectation from the vV^ Standard Model 
process, from physics processes other than uV^ , the expected contamination from cosmic rays and 

beam-related backgrounds and the overall number of expected events in the vV^ channel. The uncer- 
tainties reflect either the Monte Carlo statistics or the statistics of the control samples used for the 
instrumental background estimate. 
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Systematic errors on the cross section 


qq7 channel 


Aa (pb) 


vv^ channel 


Aa (pb) 


Efficiency 


0.115 


Efficiency 


0.045 


Jet Modelhng and Reconstruction 


0.074 






Background 


0.066 


Background 


0.006 


Cut 


0.018 


Cut 


0.005 


Energy Scale 


0.009 


Energy Scale 


0.003 


Energy Resolution 


0.008 


Energy Resolution 


0.012 


Beam Energy 


0.003 


Beam Energy 


0.001 


Luminosity 


0.020 


Luminosity 


0.005 


Total Systematic Error 


0.154 


Total Systematic Error 


0.048 



Table 4: Systematic errors on the measurement of the cross-section for the qq7 and vv^ channels. 



NTGC Fit Results and Expected Sensitivity 




qqj channel 


vu^ channel 


Coupling 


Fit Results Sensitivity 


Fit Results Sensitivity 


hf 
% 


0.000 ±0.122 ±0.26 
0.000 ±0.081 ±0.17 
-0.055lg:i^^ -0.33, +0.19 
0.035l|J:|j^i -0.14, +0.21 


0.000 + 0.166 +0.35 
0.000 + 0.115 +0.25 
-0.107 + 0.169 -0.48, +0.37 
0.067 + 0.116 -0.19, +0.25 


K 

"-2 


0.000 + 0.074 +0.15 
0.000 + 0.049 +0.10 
-0.06llg:|j|| -0.033, +0.031 
0.049lg:g^^ -0.023, +0.025 


0.000 + 0.100 +0.21 
0.000 + 0.069 +0.15 
-0.163t|j;?|5 -0.068, +0.064 
O-lSStol^l -0.045, +0.050 



Table 5: Best fit values of the anomalous couplings and the corresponding statistical errors for 
the qq7 and uVj channels. The expected statistical errors, obtained under the hypothesis of perfect 
agreement with the Standard Model expectations, are also shown. 
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Combined Results on NTGC 


Coupling 


Cc'ulral Vtihic 


95% C.L. 


Jul ('r\'al 




O.OOOiO.lOO 


[ -0.190, 


+ 0.190 ] 




0.000±0.068 


[ -0.128, 


+ 0.128 ] 






[ -0.269, 


+ 0.119 ] 




U.U^D_QQgg 


[ -0.084, 


+ 0.175 ] 




0.000±0.061 


[ -0.115, 


+ 0.115 ] 


hi 


0.000±0.041 


[ -0.077, 


+ 0.077 ] 


hi 


-0.08010-039 


[ -0.164, 


- 0.006 ] 


hi 


0.064t0;0i 


[ +0.007, 


+ 0.134 ] 



Table 6: Central values and 95% Confidence Level intervals for the eight anomalous couplings as 
determined from the combination of the qq7 and the vu^ channels. Systematic uncertainties have 
been incorporated in the errors and in the 95% C.L. limits. 



Systematic errors on NTGC 




qq7 channel 


uu^ channel 


Ahf 


0.015 


0.021 


Ahf 


0.009 


0.015 


Ahf 


0.023 


0.025 


Ah| 


0.013 


0.021 


Ah7 


0.009 


0.012 


Ah] 


0.006 


0.009 


Ahl 


0.021 


0.050 


Ah] 


0.016 


0.044 



Tabic 7: Absolute systematic errors on the NTGC couplings determined in the qq7 and in the 
fV^ channel. 
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Figure 1: (a) Distribution of the angle between the photon candidate and the closest track; the sharp 
edge at 0.26 rad reflects the size of the isolation cone, (b) Distribution of the total charged energy 
in the isolation cone; the explicit cut on this variable is set to 2 GeV. In both cases, the distribution 
refers to qq7 events which are accepted or which fail only a single selection criterion. The Monte 
Carlo prediction is normalised to the integrated luminosity of the data. 
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Figure 2: Distribution of E^, cos a-y-jct and cos 9-y measured in qq7 events. The data (dots) are super- 
imposed on the Standard Model predictions from the KK2f Monte Carlo (solid histogram), normalised 
to the integrated luminosity of the data sample. The shaded area in the histograms represents the 
background from four-fermion production and two-photon interactions and the contamination from 
photon candidates not corresponding to truly initial state radiation in e'''e~— >qq events. 
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Figure 3: Distribution of kinematic variables and cos^^ measured in vV^ events. The data (dots) 
are superimposed on the Standard Model prediction (solid histogram), taken as the average of the 
two generators KORALZ and NuNUGPV98, normalised to the integrated luminosity of the data sample. 
The negligible (<0.3%) background from physics processes other than vV^ is not shown. 
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Figure 4: Distribution of kinematic variables E^, cosa^_jet and cos 9^ measured in qq7 events in the 
data (dots) and in the KK2f Monte Carlo (solid histogram). The simulated Standard Model expecta- 
tions, reweighted to incorporate the contributions from hg = ±0.5 are also shown. The normalisation 
of the Monte Carlo predictions is defined by the number of events selected in the data. 
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Figure 5: Distribution of kinematic variables and cos 9^ measured in 1^17^ events in the data (dots), 
and in the Standard Model Monte Carlo (solid line). The simulated Standard Model expectations, 
reweighted to incorporate the contributions from hg = ±1 are also shown. The normalisation of the 
Monte Carlo predictions is defined by the number of events selected in the data. 
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Figure 6: Negative log-likelihood function for the hf couplings as obtained from the analysis of the 
qq7 channel (dash-dotted line), of the lyi'j channel (dashed line) and from their combination (solid 
line) . 
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Figure 7: Negative log-likelihood functions for the hj couplings as obtained from the analysis of the 
qq7 channel (dash-dotted line), of the lyi'j channel (dashed line) and from their combination (solid 
line) . 
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